The Average Kinetic Energy of the Heavy Quark in A5 
in the Bethe-Salpeter Equation Approach 



X.-H. Guo"i and H.-K. Wu^^ 
. Institute of Low Energy Nuclear Physics, Beijing Normal University, Beijing 100875, China 
b. Physics Department, Louisiana State University, Baton Rouge, LA 70803, USA 

Abstract 

In the previous paper, based on the SU{2)f x SU{2)s heavy quark symmetries of 
the QCD Lagrangian in the heavy quark limit, the Bethe-Salpeter equation for the 
heavy baryon A;, was established with the picture that A^ is composed of a heavy 
quark and a scalar light diquark. In the present work, we apply this model to calcu- 
late /i^ for Aft, the average kinetic energy of the heavy quark inside A^. This quantity 
is particularly interesting since it can be measured in experiments and since it con- 
tributes to the inclusive semileptonic decays of A;, when contributions from higher 
order terms in l/M^ expansions are taken into account and consequently influences 
the determination of the Cabibbo-Kobayashi-Maskawa matrix elements Vuh and Vcb- 
We find that /i^ for A;, is 0.25GeV^ ~ 0.95GeV^, depending on the parameters in 
the model including the light diquark mass and the interaction strength between the 
heavy quark and the light diquark in the kernel of the BS equation. We also find 
that this result is consistent with the value of fi^. for A;, which is derived from the 
experimental value of /i^ for the B meson with the aid of the heavy quark effective 
theory. 
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I. Introduction 



The physics of heavy quarks has attracted intense interests in recent years, partly 
because of the discovery of the flavor and spin symmetries in QCD, SU (2) / x SU (2)s, 
in the heavy quark limit and the establishment of the heavy quark effective theory 
(HQET) [1]. Compared with the research on heavy mesons, heavy baryons have 
been studied less both theoretically and experimentally. However, more and more 
experimental data for heavy baryons have been and will be obtained. This will help 
to test theoretical predictions for heavy baryons. For example, the lifetime of A;, 
has been measured in several experiments [2] . The measurement of the nonleptonic 
decay of A;,, A;, AJ/'^, has been done [3]. There have also been the measurements 
of the semileptonic decays of A^, A^ — > Ac/ anything [4] and A;, AJ~i/i [5]. On 
the other hand, since heavy baryons are composed of three quarks instead of two, 
theoretical studies for heavy baryons become more complicated. In order to under- 
stand the hadronic structure of heavy baryons, more theoretical and experimental 
studies are needed. 

HQET can simplify the physical processes involving heavy quarks. A^ is com- 
posed of a heavy b quark and two light quarks, u and d. When the heavy quark 
mass is very large compared with the QCD scale Aqcd, the light degrees of freedom 
(the light quark system) in a heavy baryon becomes blind to the flavor and spin 
quantum numbers of the heavy quark because of the SU{2)f x SU{2)s symmetries. 
Therefore, the angular momentum and flavor quantum numbers of the light degrees 
of freedom become good quantum numbers. Hence it is natural to regard the heavy 
baryon A;, to be composed of a heavy quark and a light scalar diquark, [ud]o, with 
[ud] flavor quantum number and zero spin and isospin. 

The Bethe-Salpeter (BS) equation is a formally exact equation to describe the 
relativistic bound state [6, 7, 8]. In the heavy quark limit the BS equation can be 
simplified to a great extent and has been applied to give many theoretical results 
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concerning heavy mesons and heavy baryons [9, 10, 11, 12, 13, 14]. With the model 
for the composition of which is described above the heavy baryon A;, is reduced 
from a three-body system to a two-body system. In this picture the BS equation 
for Aft was estabhshed [10, 11]. A scalar confinement and a one gluon exchange 
term compose the kernel of the BS equation in this model. Furthermore, this model 
was generahzed to the heavy baryons E^*^ , S^*^ , and Q^*^ which are regarded to be 
composed of a heavy quark and an axial- vector diquark [12]. 

In HQET the strong interaction of a heavy quark with four-velocity v can be 
described by the following Lagrangian density [1, 15]: 

jC = Kiv ■ DK + 7^K[{iD^f]K + -^Ka^uG^^K, (1) 

ZMq ^Mq 

where denotes the field of the heavy quark, Mq is the mass of the heavy quark, 
^ Qn _ ig^A^ is the covariant derivative, D± — — v^v ■ D, and G^^ is the 
gluon field tensor. The second operator in Eq. (1) is related to the average kinetic 
energy of the heavy quark due to the residual motion of the heavy quark inside the 
heavy hadron and the third one corresponds to the spin energy of the heavy quark. 
The kinetic energy and the spin energy of the heavy quark can be described by the 
following two local matrix elements respectively: 

2 _ {HQ\K{iDA.fK\HQ) 
/^.- 2M ' 

and 

2 _ {HQlgghya^^Gfj^iyhy \ Hq) 
^« " AM ' 

where Hq {Q — h or c) denotes a heavy baryon containing a heavy quark Q and M 
is the mass of Hq. 

The parameters /x^ and /x^ are of particular interests since they contribute to 
the inclusive semileptonic decays of heavy hadrons when contributions from higher 
order terms in 1/Mq expansions are taken into account and, therefore, influence the 
determination of the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements Vub and 



Vcb- Therefore, it is very interesting to calculate these nonperturbative quantities 
theoretically. 

There have been extensive studies in literature on inclusive semileptonic decays 
of bottom hadrons, Hh XeUe, especially since the establishment of HQET [16, 
17, 18, 19, 20, 21, 22]. These studies include corrections to the leading order results 
both from perturbative QCD (^^(Mi,)) terms and from nonperturbative terms which 
are suppressed by powers of M^. It has been pointed out that there is no I /Mb 
corrections to the leading order result in l/M^ for the differential decay width of 
semileptonic decays of bottom hadrons, dV/dq^dEg, where q is total momentum 
of the electron and the neutrino and Eg is the electron energy [16]. Then Bigi et 
al. studied 1/M^ corrections to the decay width dT/dEg [21]. Manohar and Wise 
analyzed extensively corrections to dV/dq^dEe for unpolarized bottom hadron 

Hb and for polarized A;, [22]. In recent years, theoretical calculations for the inclusive 
semileptonic decay widths and for the moments of inclusive observables have been 
carried out to order l/M^ and q;^/3o (/3o = H - 2n//3) [23, 24, 25, 26, 27]. It was 
found that the 1/M^ corrections are characterized by the two parameters /x^ and 
fiQ, which can be extracted from experimental data and theoretically should been 
determined in a nonperturbative way. 

In the case of the B meson, the parameter /x^ can be extracted from the data for 
the hyperfine splitting between B and B* mesons. The other parameter /x^ has been 
extracted from the experimental data for the inclusive semileptonic B meson decays 
B XJP and B Xs^f [28, 29, 30, 31, 32, 33]. In Ref. [33], using the theoretical 
formulae provided in Refs. [24, 25] , Buchmiiller and Flacher obtained the most recent 
result for /x^, //^ = 0.401 ± O.OAOGeV"^, from a combined fit to the moments of the 
hadronic mass distribution and the moments of the leptonic energy spectrum in 
B — > XJi/ and the moments of the photon energy spectrum in 5 — > Xg'j which 
are measured in the BABAR, Belle, CDF, CLEO, and DELPHI experiments [32, 
34, 35, 36, 37]. Theoretically, /x^ has been calculated in various phenomenological 
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models such as QCD sum rules and the BS equation [14, 38, 39, 40, 41] and by lattice 
QCD [42]. The theoretical results of /^^ for the B meson depend on models strongly. 
Some of them are consistent with the experimental value, //^ = 0.401 ± 0.0406*6^^, 
while some of them, including that from the BS approach for the B meson, are not. 
This needs further and more careful investigations. 

Compared with the case of the B meson. A;, has been studied less both exper- 
imentally and theoretically. Since is composed of a heavy quark and a light 
scalar diquark, the parameter is zero for A^. Although there has been no direct 
experimental measurement of /i^ for A^, one can expect it to be measured in the 
future since more and more data on will be collected. Furthermore, with the 
aid of HQET, /^^ for A^ can be related to /^^ for the B meson [22]. Hence, one can 
derive the value of /i^ for A;, from the experimental value of //^ for the B meson. 
Therefore, it is important to give results for //^ for A^ from theoretical calculations. 
The aim of the present work is to calculate the average kinetic energy of the b quark 
in the heavy baryon A^ with the BS equation model for A^ [10, 11]. We will give 
the numerical result for this parameter, discuss its dependence on the parameters 
in the model, and compare our result with the value of for A^ derived from the 
experimental value of //^ for the B meson through HQET. 

The reminder of this paper is organized as the following. In Section II we review 
the basic formalism for the BS equation for A^. In Section HI we give numerical 
solutions for the BS wave function and then apply the BS equation to calculate 
nl. numerically. We also discuss the dependence of our result on the parameters 
in the model and compare this result with the value of /^^ for A^ derived from the 
experimental value of //^ for the B meson with the aid of HQET. Finally we give a 
summary and discussion in Section IV. 

II. Formalism for the BS equation for A;, 

As discussed in Introduction, Aq is regarded as the bound state of a heavy quark 



5 



and a light diquark. Based on this picture the BS wave function of Aq is defined as 
follows: 

X(X1,X2,P) = (O|rV'(^l)0(^2)|AQ), (4) 

where ip{xi) and 4>{x2) are field operators of the heavy quark and the diquark, 
respectively, and P is the momentum of Aq. The BS wave function in the momentum 
space, Xp(p), is related to x{^ij^2, P) through the following equation: 

X{x^,x2, P) = e^^^ / ^;^p(p)e^^>-, (5) 

where p and x{= xi — X2) are the relative momentum and the relative coordinate 
of the heavy quark and the light scalar diquark, respectively, and X is the center 
of mass coordinate which is defined as X — XiXi + X2X2, where Ai = m^Md ' 
A2 = m^+Mq ' "with Md being the mass of the diquark. The momentum of the heavy 
quark is pi = AiP + p and that of the diquark is p2 = —X2P + P- 
The mass of the heavy baryon, M, satisfies the following relation: 

M^Mq + Md + Eo + 0{^), (6) 

Mq 

where Eq is the binding energy in the leading order of 1/Mq expansion. 
The BS equation in the momentum space can be written as follows [10]: 

Xp(p)^Sf(XiP + p) J ^,G{P,p,q)xp{q)Sr>{-X2P + p), (7) 

where G{P,p,q) is the kernel which is defined as the sum of all the two particle 
irreducible diagrams with respect to the heavy quark and the light diquark. Sp 
and Sd in Eq. (7) are propagators of the heavy quark and the light scalar diquark, 
respectively. 

The kernel G{P,p, q) includes two terms in the model: a scalar confinement term 
Vi and a one gluon exchange term V2 [9, 10, 12], 

-iO^I® IVi + v^(» {P2 + p'2YV2, (8) 
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where p2 and the momenta of the hght diquark attached to the gluon. The 

vertex of the gluon with the diquark depends on the structure of the diquark. This is 
taken into account by introducing a form factor -^[(^2—^2)^]) which is parameterized 
as F{Q'^) = ^'a [43], where Ql is a parameter which freezes FiQ"^) when Q'^i— 
ipt — QtY) is very smalL 

It has been shown that in the leading order of 1/Mq expansion we only need 
one scalar function, (pp^p), to describe the BS wave function [10]. (I)p{p) is related 
to Xp{p) the following: 

Xp{p) = (ppiphAgiv), (9) 

where v is the velocity of the heavy baryon and u\q (v) is the spinor of the heavy 
baryon. 

Define the longitudinal and transverse momenta with respect to v: pi — v ■ p — 
A2M, pt = p — {v ■p)v. Using the covariant instantaneous approximation, pi = qi, at 
the vertex of the heavy quark and the gluon, we have the BS equation in the leading 
order of 1/Mq expansion, 

where (j)pipt) = f{dpi/'^'^)<ppip) and V stands for V in the covariant instantaneous 
approximation pi — qi. 

Integrating Eq. (10) by / dpi/2TT and applying the residue theorem we obtain 
the equation for the BS wave function, ^p{pt), 

= 2WA-W~^E,^Mo) / - '^.2>4)^P(„). (11) 

The kernel Vi and V2 have the following expression in the case of the heavy 
baryon [10, 12]: 



where k and af^ are coupling parameters related to the scalar confinement and 
the one gluon exchange diagram, respectively. The second term in Eq. (12) is the 
counter term which removes the infra-red divergence in the integral equation. The 
parameter u is introduced to avoid the infra-red divergence in numerical calculations. 
The limit m — > is taken in the end. 

Substituting Vi and V2 into Eq. (11) we have 



2Wp [J 47r2 (p| + g2 + ^2)2 _ 4^2,^2 

/ 



^ 327r{affQlW, f q^dq^ 1 



3{Qo-u^) J 47r2 2ptqt 
2 I /n2 , 

0p(% 



{pt-qtY + u^ 



_ (Pt + qt? + Ql 



ipt-qtY + Ql 



qjdqt 
47r2 



(Pt + qt + uv - Mqi 



III. Calculation of the average kinetic energy of the b quark 

in Ab 

In this section we solve the BS equation numerically and then apply the results 
to calculate the average kinetic energy of the b quark inside the heavy baryon A^,, 
nl., which is defined in Eq. (2). The BS wave function for A;,, (ppipt) in Eq. (14), 
can be solved numerically by discretizing the integration region (0, 00) into n pieces 
{n is chosen to be sufficiently large). We use the n-point Gauss quadrature rule to 
evaluate the integral. ThenEq. (14) becomes an eigenvalue equation. The numerical 
results for 0p(pt) are obtained by solving this eigenvalue equation. Eq. (14) is 
a homogeneous equation which leaves the normalization of 0p(pt) undetermined. 
We use the following normalization condition to fix the amplitude of the BS wave 
function^: 

(AblJi^/i^lAb) = 2M. (15) 



^One can also use the expression = _ to calculate /x^. This expression is 

independent of how the BS wave function is normalized. 
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In the model we have several parameters, i.e. a'f, k,, Ql, Md and Eq. The param- 
eter Ql is taken as Ql ~ 3.2GeV^ [10, 12, 43]. The parameters a'f and k are related 
to each other when we solve the eigenvalue equation with a fixed eigenvalue [10]. 
The parameter k varies in the region between 0.02GeV^ and O.lGeV^ [10]. From Eq. 
(6) the parameters and Eq are constrained by the relation Mo + Eq = M — 
for Aft in the leading order of l/M^ expansion. In our numerical calculations we use 
Mb = 5.02GeV which leads to consistent predictions with experiments from the BS 
equation in the meson case [9]. Consequently we have Md + Eq = 0.62GeV for A;, 
(where we have neglected 1/Mi, corrections). The parameter can not be deter- 
mined and hence we let it vary within some reasonable range. For A^, we choose 
Mo to be in the range 0.65GeV ~ O.SOGe^. With this choice for Md, the binding 
energy Eq is negative and varies from around —SOMeV to —ISOMeV. 

The numerical results for af^ corresponding to various values of k are given in 
Tables 1, 2, and 3. Then the numerical results for the BS wave function depend 
on two parameters, k, and Md- In Fig. 1 we show the solutions for the BS wave 
function for some typical values of k and Md- 

Since /x^ is a Lorentz scalar [22] we are free to choose a special frame for the 
calculation of this parameter. For simplicity we choose the rest frame of A^ in which 
Eq. (2) becomes 



.2 _ {AQ\h,:{iDfh,]AQ) 



2M ^ ^ 

The diagram for calculating the average kinetic energy of the b quark inside A^ 

is shown in Fig. 2. Assuming the light diquark acts as a spectator, we obtain the 

following expression for //^ which is related to the BS wave function of A^: 

A = l^j -^,Xp{p)fxp{p)SD\-\2P + p). (17) 

Substituting Eq. (9) and the relation between 0p(p) and ^p{pt), Eq. (10), into 
Eq. (17) and integrating the pi component by selecting the proper contour we have 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

(GeV) 

Fig. 1 




The three-dimensional integral in Eq. (18) can be simplified to one-dimensional 
integral. This leads to 

^'-2¥ (^^) 

As shown in Fig. 1, the numerical results for the BS wave function (j)p{pt) depend 
on the parameters k (or a^^) and Mjj. Therefore, the results for /i^ also depend on 
these parameters. For example, taking Mid = 0.7GeV and k = O.OAGeV^, we get 
/i^ = OATGeV"^. In Tables 1, 2, and 3 we list the numerical results for /i^ for various 
values of the parameters k, and Mjy. 
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It can be seen from these tables that the value of /x^ changes from 0.25GeV'^ 
to 0.95GeV^ in the variation ranges of the model parameters Md and k. The 
dependence of the average kinetic energy of the b quark inside on the model 
parameters is quite strong. Furthermore, the dependence on k is stronger than 
that on Ml). For instance, for M/j = 0.7GeV, when k varies from 0.02GeV^ to 
O.lGeV^, the change of is about OAAGeV"^; for k, = OMGeV^, when Md varies 
from O.GbGeV to O.SOGeV, the change of fil is about 0.27GeV'^. 

In principle, the parameters in the model can be determined through the com- 
parison between theoretical predictions and experimental measurements about some 
physical processes of if the data are accurate enough. In Refs. [10, 11, 13] some 
phenomenological predictions for A^ such as semileptonic and nonleptonic decay 
widths of Aft are given in the BS equation approach. Since the heavy quark mass 
is not infinite in reality, in order to give more exact phenomenological predictions 
1/Mq corrections to the BS equation for Aq are analyzed in Ref. [11] based on 
the assumption that Aq is composed of a heavy quark and a scalar light diquark. 
Including both the 1/Mq corrections and the QCD corrections to the weak decay 
form factors [44] the prediction for the decay rate for A& AcW can be obtained 
as 2.70 ~ 4.07 x 10^°s~^ in the variation ranges of Md and k, (where Vcb is taken as 
0.042 [45]). The uncertainty of this prediction is mostly from the uncertainty from 
K. The experimental data for the decay rate for A;, — > AJp, which is in the range 
2.3 ~ 6.7 X 10^°s~^ [5], is consistent with the prediction. Therefore, we can not 
determine the parameters in the BS equation model from this process at present 
due to the large error in the experimental data. With more and more data avail- 
able in the future one can constrain the parameters in the BS model much better. 
Furthermore, the experimental data for the nonleptonic decay widths for A;, — > Ac 
plus a pseudoscalar or a vector meson (the predictions for them have been given in 
Ref. [11]) can also be used to determine the parameters in the model. 

Theoretically, there have been some phenomenological calculations for the di- 
ll 



quark mass from the BS equation for the diquark [46] and from the relativistic 
potential model for the diquark [47], respectively. The masses for the [ud]o diquark 
obtained in these two approaches depend on the model parameters and are consistent 
with what are used in our BS model. 

As mentioned in Introduction, although there has been no direct experimental 
measurement of /i^. for A^, one can relate this quantity to /i^. for the B meson with 
the aid of HQET [22] . In this way, one can derive the value of //^ for Af, from the 
experimental value of for the B meson. It was shown that when the masses of 
heavy hadrons are expanded to order 1/Mq one has the following relation: 

f^'M-f^UB) = ^^§^^^^{[M(A,)-M(i^),g-[M(A,)-M(S),g}, (20) 

where i^l(Ab) (i^UB)) is i^l for A^ (B), M{B) (M{D)) is the mass of B (D), and 
M(D)avg (M(i?)avg) is defined as the spin averaged mass of D (B) mesons (for 
instance, M{D)^yg = [M{D) +3M{D*)]/A for D mesons). The masses of D, B, and 
Ac have been measured accurately and the largest uncertainty of the right hand side 
of Eq. (20) comes from the mass of A;, [45]. Using M{Ab) = 5624 ± 9MeV and the 
masses of D, B, and Ac provided in Ref. [45], the right hand side of Eq. (20) is 
0.025 ± 0.052Gey^ where the error comes mostly from the error of the mass of A^ 
(the errors of the masses oi D, B, and Ac contribute httle). Consequently we obtain 
uliK) from uliB) = 0.401 ±0.040Gel^^ (which was obtained by fitting the data in 
the so-called kinetic scheme [33]) as follows: 

= 0.426 ± 0.066G'ey^ (21) 

where the error includes those from both //^(-B) and the mass of A^. 

Besides the uncertainty in Eq. (21), the 1/Mq terms in the expansion for the 
masses of heavy hadrons may also cause some uncertainty to /^^(A;,). Two param- 
eters, p\) and p^, appear in the 1/Mq terms in the masses of A^, Ac, and the spin 
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averaged masses of D and B mesons'^ [40, 48]. The parameter p|, has been ex- 
tracted from the fit in Ref. [33] while p^, which is a nonlocal correlator of the two 
operators h-t,{a ■ -D)^/i^, has not been determined, is of order (A is defined 
as the difference between the mass of a heavy hadron and the mass of the heavy 
quark inside the hadron in the heavy quark limit) [33]. p^ is also expected to be 
of order A^. Although there may be some cancellation between the parameters p^ 
and p^ for the heavy baryons and those for the heavy mesons in the mass difference 
[M(Ae) - M(L))avg] - [M(Ab) - M(5)avg] on the right hand side of Eq. (20), we 
assume that the 1/Mq terms in this mass difference is of order A^/Mq to make a 
conservative estimate on the influence of the 1/Mq terms on /i^(Af,). The h? /M^ 
terms give the main contribution to A*^(A(,) in Eq. (20), which is about Q.Q^GeV^ if 
we take A to be Q.QGeV [33]. 

Taking into account all the uncertainties from p^B), the mass of A^, and the 
A? /M^ terms in the masses of heavy hadrons, one may expect A*^(Ab) to be roughly in 
the range Q.21GeV'^ ~ O.SSGeF^. This is consistent with our result in the BS model, 
0.25GeV'^ ~ 0.95GeV'^. Conversely, one may give a rough constraint on the ranges 
of the parameters in the BS model from the range of /x^(Ab), 0.27GeV^ ~ 0.58GeV'^. 
For instance, when Md is 0.65GeV, k is roughly in the range 0.02GeV^ ~ O.OSGeV^ 
from Table 1, while when Md are 0.7GeV and O.SGeV, k are roughly in the ranges 
0.02^6^3 - O.OGGeV^ and Q.02GeV^ - QMGeV^ from Tables 2 and 3, respectively 

IV. Summary and Discussion 

The average kinetic energy of the b quark inside A^, p"^, is an interesting quantity 
both theoretically and experimentally. It contributes to the inclusive semileptonic 

decays of A;, when contributions from higher order terms in 1/Mh expansions are 

taken into account and influences the determination of the CKM matrix elements 

contain four terms, p^^, p^Q, p|, p\, while only p^^ and p| contribute to the masses of A;, 
and Ac and the spin averaged masses of D and B mesons [40, 48]. 
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Vub and Vcb- By comparing the experimental data with the theoretical predictions 
for such decays one can extract the value of /i^. 

Based on the BS equation model for the heavy baryon Ab, which is regarded as 
composed of the heavy b quark and a light diquark, we have calculated the average 
kinetic energy of the b quark inside A^. The kernel of the BS equation consists of a 
one gluon exchange term and a scalar confinement term. Since /x^ is expressed as 
the overlap integral of the BS wave function of A^, we first solved out this BS wave 
function numerically by transfering the integral equation for the BS wave function 
into an eigenvalue equation. We have found that the value of /i^ varies in the region 
between 0.25GeV'^ and 0.95GeV^ depending on the parameters in the model. The 
dependence of /x^ on the parameters in the model was discussed in some detail. 
We have compared our result with the value of /i^ for A^ which is derived from the 
experimental value of fi^ for the B meson with the aid of HQET and found that they 
are consistent. Conversely, the latter may also be used to give a rough constraint 
on the parameters in the BS model. 

Compared with the meson case, heavy baryons are much more complicated since 
there are three quarks in a baryon. Even though we have simplified the bound 
state equation for a heavy baryon with the diquark picture, large uncertainties are 
still introduced in the BS equation for the heavy baryon. This is reflected in the 
large ranges of the parameters in the model, i.e. k. and Mjj. This leads to a much 
larger range for the phenomenological prediction for the average kinetic energy of 
the b quark inside A^. Fortunately much more data will be available in the future 
experiments, e.g. LHCb. This provides an opportunity to constrain the model 
parameters more accurately by comparing the experimental data with the BS model 
predictions for the physical processes, say semileptonic and nonleptonic decays of 
Afc. 
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Figure captions 

Fig. 1 Numerical results for the BS wave function (j)p{pt)- The solid (dashed) line 
corresponds to Md = Q.lGeV and n = 0.02(0.1)6*6^^ The dotted (dot-dashed) 
line corresponds to k = ^MGeV^ and Md = 0.65(0.8)G'ey. 

Fig. 2 The diagram for calculating the average kinetic energy of the b quark inside 
Afo. The black dot represents the operator hy{iDYhy. 





0.02 


0.04 


0.06 


0.08 


0.10 


af 


0.62 


0.67 


0.70 


0.72 


0.75 




0.25 


0.39 


0.51 


0.62 


0.72 



Table 1: The values of af^, and the corresponding //^ for Md — 0.656*6^. 
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Table 2: The values of k, af, and the corresponding /^^ for Md = 0.7GeV. 
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Table 3: The values of af, and the corresponding //^ for Md — O.SGeV. 
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